Rise or sink: spherical intruder in density 
frequency dependent static granular fluid 
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Abstract 

A simple model for intruder in vibrating granular bed is con- 
structed by introducing a term that governs frequency dependent 
granular bed density. Varying the vibrating frequency will drive the 
buoyant force acting on the intruder by the granular bed. Swapped 
regions of rising and sinking intruder compared to other reported re- 
sult have bee found. 

Keywords: frequency dependence, granular bed, spherical intruder. 



1 Introduction 



A phenomenon known as Brazil nut effect (BNE) is introduced descriptively 
by a work via molecular dynamics (MD) simulation for one and several in- 
truders PQ. The density difference between intruders and granular bed plays 
important role in size separation of the grains [2] , which determines whether 
a larger intruder could sink instead of rising [3J. Vibrating amplitude, for 
binary mixture, could drive the state of the mixture into BNE or its reverse 
(reverse Brazil nut effect, RBNE) [H [5]. It is interesting to consider that 
granular bed acts as fluid when the vibrating amplitude and frequency are 
not zero, but when these parameters are zero, it should be a solid. A simple 
model that treats the granular bed as fluid is presented with modification 
to the granular fluid density, which is previously independent from vibrat- 
ing amplitude and frequency [6]. This model is simpler than that is already 
reported [7J. 



2 Theory 

A spherical intruder with radius R and density p is placed in a granular fluid 
with density p g . As it moves updward or downward in the granular fluid it 
will be under influence of earth gravitational force Fq 

f g = -\*pgR\ (i) 

buoyant force F B 

F B = ^ Pg gR\ (2) 

and viscous drag force Fp 

F D = -Q*vR^, (3) 

with upward y-direction is taken to be positive. Constants g and rj repre- 
sent earth gravitational acceleration and granular fluid viscosity, respectively. 
There are no additional complex forces, such as thermal buoyancy force, com- 
presive force, and dynamic tensile force [TJ. 
Granular fluid density is defined as 
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Pg = (p + Por)exp(-r), (4) 

where p is density of one particle of granular bed and T is dimensionless 
acceleration 

r = W (5) 
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with A is vibration amplitude and / is vibration frequency, which is common 
used in indicating the vibration influence [H EJ [7]. 

Using Newton's second law of motion with Equation ([I])-© will give a 
differential equation 

, d y , n fet\ 



with 
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47r[p-(p + p r)exp(-r)]fi 3 

C2 " 3^ • (8) 



Equation (jf)]) has solution 
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where Vq and yp are initial velocity and position, respectively. Since the 
intruder always has no initial condition or vq = 0, then Equation (jUJ) will be 
reduced to 

y(t) = (|) exp (- Cl t) +(|)*+ U - |) • (10) 

3 Results and discussion 

Because of physical properties of intruder and granular bed Equation ([Zj) 
restricts that c% > 0. For C2 it can be positive or negative according to 
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Equation (jSJ) since the intruder can have larger or smaller density value than 
density of one particle of granular bed. 

From Equation f JTOj) velocity of the intruder can be found, which is 

«(*)= (|) [1- exp (-erf)], (11) 

which tells us that there is a terminal velocity that can have positive and 
negative value. Positive and negative value means a BNE and a RBNE. 
respectively. Using Equation (J7J) and (jSJ), Equation ffTTl) can be rewritten 
in such form that shows the influence of density difference between intruder 
and one particle of granular bed, which is 



W (t)oc[p-(p + p r)exp(-r)], (12) 

which shows that whether the intruder will rise or sink is dependent ex- 
plicitely on T. 

Figure [1] shows the room parameter of d/do againts m/mo where the re- 
gions for BNE and RBNE are swapped as reported in [5] , but the dependence 
of crossover line to T is similar to what is found by other [5] . As F increases 
the gradient of transition line in room d/do againts m/rriQ also increases. For 
T about 8 in the range used in Figure [1] gives no RBNE states (not shown). 

Rise time for the intruder can be found by setting initial and final value 
of y in the Equation (flOl) . In this case y(0) = yo and y{r) = yf will be used 
to obtain the rise time r. 



y/= (|) exp( - cir)+ (l) r+ r-| 

^y/-yo + |=(|)exp(- Cl r) + (|)r 
=> 1 + f^J (Vf ~ Vo) = cit + exp (-cir). 



(13) 



Considered that there is r where r >> r will reduce Equation ffl3|) into 

1 , fci 



-+ - (Vf-Vo), (14) 

Ci VC 2 , 



which can be written in 
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r+ ~ ex + [ P - ( P + Po r) ex P (-r)] ' (15) 

with 

_ 3m Cl (y f -y ) 

C3 " 4^ ' (16) 

And for r << r will reduce Equation ( Fl3|) into 

r _~_Il n (l + \ (17) 

ci I [p-(p + p r)exp(-r)]i 1 ; 

Figure [2] shows unexplained results that there are negative rise time for 
both calculations using Equation ( TT5|) and (|T7|) . There are two things to do, 
solve Equation ( fT3l) directly without approximation or redefine Equation (j4j) 
for the granular bed density. 



4 Conclusion 

A model treating granular bed as static fluid has already done. Room pa- 
rameters for BNE and RBNE shows swapped regions as reported by other. 
Rise time has also negative values. Further investigation is needed. 
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Figure 1: Room parameter of d/d againts m/m for T: (a) 0.5, (b) 1, (c) 2, 
and (d) 4. BNE and RBNE are indicated by o and □, respectively. 
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Figure 2: Rise time for t_ (o) and r + (□). 



